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Abstract. The rare B decays B° (s) -> B° -> K* /u + fi~ and B° s -> 0y are studied using up to ~ 0.41 fb -1 of 

pp collisions at = 7 TeV collected by the LHCb experiment in 2010 and 201 1. A search for the decays 5?x — > 

is performed with 0.41 fb -1 . The absence of significant signal leads to &(B° S — > < 1.4 x 10~ 8 and 

£>(B° — > jj, + /jT) < 3.2 x 10~ 9 at 95 % confidence level. The forward-backward asymmetry, fraction of longitudinal 
polarization and differential branching fraction of B° — > K*°jj, + /j,~ , as a function of dimuon invariant mass, are 
measured in 0.31 fb -1 . The ratio of branching ratios of the radiative B decays B° — > K*°y and B? — > (py has been 
measured using 0.34 fb" 1 . The obtained value for the ratio is 1.52 + 0.14(stat) + O.lO(syst) + 0A2(f s /f d ). Using 
the HFAG value for &(B° -> K*°y), !B(B° S -> <fry) has been found to be (2.8 + 0.5) x 10" 5 . 



1 Introduction 

The LHCb experiment |1 ] has provided preliminary re- 
sults in the measurement of the forward-backward asym- 
metry, fraction of longitudinal polarization and differential 
branching fraction of B° — > K*°jj + fi~ Q and the measure- 
ment of the &(B° S -> 0y) 0. LHCb has also provided up- 
per limits in 8(BP S -> /z + /T) and !B(B° -> //+//") |4|. Sect. [2] 
sumarizes the analysis and results obtained by LHCb in the 
study of B° — > K*°/i + /d~ . Sect, p] sumarizes the measure- 
ment of !B(B° S ^ (f>y)/!B(B ^ K^y) and Sect. [^sumarizes 
the analysis and results of #9 , — > . 



2£ ^ K*°fi + fi- 

The rare decay 5° — > K*°/j, + /u~ is a Z? — > s, flavour chang- 
ing neutral current decay, mediated by electroweak box 
and penguin diagrams in the Standard Model (SM). In mod- 
els beyond the SM, new particles can enter in competing 
loop-order diagrams resulting in large deviations from SM 
predictions (see for example Refs. 15161 ). 

B° — > K*®ia + ia~ candidates are selected by first apply- 
ing a loose pre-selection based on the B° lifetime, daugh- 
ter impact parameters and a requirement that the B° points 
back to one of the primary vertices in the event. A tighter 
multivariate selection, based on a boosted decision tree 
(BDT), is then applied to select a clean sample of B° — > 
K*°ju + ji~ candidates, with a signal-to-background ratio in 
a 100MeV/c 2 window around the reconstructed B° mass 
of about three-to-one. The BDT is based on the B° kine- 
matics, B° vertex quality, daughter track quality, impact 
parameter and kaon, pion and muon particle identification. 
The offline selection criteria are explicitly chosen to min- 
imise angular acceptance effects. The multivariate selec- 
tion was trained using B° d — > J/ifsK*° candidates from the 
2010 data as a proxy for the signal and B° —> 
candidates from the upper mass sideband of the 2010 data 
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for the background. Specific vetoes are used in order to 
eliminate non combinatorial background. 

The trigger, reconstruction and offline selection can all 
bias the measured angular distribution of B° — > K*°fi + fi~ 
candidates. The detection acceptance is accounted for by 
weighting events when fitting for A FB , F L and dBF/dq 2 
(where q 2 is the di-muon mass squared). Event weights are 
calculated on a per-event basis in a small phase space win- 
dow around each candidate, using fully simulated Monte 
Carlo (MC) simulation events. Simulated events are re- 
weighted to account for known data-MC differences in PID 
performance, impact parameter resolution, tracking effi- 
ciency and track multiplicity. 

The fit results for A FB , F L and dBF/dq 2 , and their com- 
parison with theoretical predictions |7 ], are shown in Fig.[T] 

The systematic error on A FB , F L and dBF/dq 2 is typ- 
ically ~ 30% of the statistical error. In the high-g 2 re- 
gion, the dominant contribution to the systematic uncer- 
tainty comes from the overall uncertainty on the accep- 
tance correction which is dictated by the limited simulation 
statistics. This can clearly be improved for future analyses. 
Throughout, a sub-dominant contribution comes from the 
data-derived performance corrections. In particular, from 
knowledge of the PID performance and tracking efficiency 
in data. This is again statistically limited and can also be 
improved with larger datasets. When fitting for A FB and 
F L the signal and background mass model and the angular 
model for the background have been varied and yield cor- 
rections at the level of 10-20% of the statistical uncertainty. 
The uncertainty on the differential branching fraction in- 
cludes the ~ 4% uncertainty coming from the measured 
B° —> J/if/K*° and J/if/ —> branching fractions 0. 

These measurements are current world best, and don't con- 
firm previous hints of a non-SM value of A FB at low q 2 . 

In the SM, the amplitude of these sy penguin transitions 
is dominated by a virtual intermediate top quark coupling 
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to a W boson. Extensions of the SM predict new heavy 
particles that may propagate virtually within the loop and 
modify the dynamics of the transition. Therefore, these ra- 
diative modes are promising laboratories that could reveal 
the presence of new phenomena beyond the SM with the 
precise measurement of the branching ratios, asymmetries 
or angular distributions. The offline selection of both the 
B° — > K*°y and B° — > (py decays is performed with 
the strategy of maximizing the cancellation of systematic 
uncertainties when performing the ratio. The analysis of 
~ 341 pb" 1 of LHCb data gives: 



&(B ^ K*°y) 
®(B° s ^<f>y) 



= 1.52±0.14(stat)±0.10(syst)±0.12(/ s // rf ) 

(1) 



Where fd (f s ) are the probabilities of the b quark to hadronize 
into B° (B° s ). This results is compatible within 1.6 standard 
deviations with the theory prediction. 

Combining the ratio of branching fractions in [T] with 
the World Average measurement for the !B(B° — > K*°y ) 
from [9], we obtain, 



S(fl?->0y) = (2.8±O.5)xl(r 



(2) 



which agrees within 1.6 standard deviations with the pre- 
vious experimental measuremen, and wich correspond to 
the most precise measurement of this BR to date. 



4 B° 



The SM predictions for the branching fractions of the FCNC 
decays B° s p + p~ and B° p + p~ are S(B° S p + p~) = 
(3.2 ± 0.2) x 10" 9 and S(B° p + p~) = (0.10 ± 0.01) x 
10" 9 ifTOl . However, contributions from new processes or 
new heavy particles can significantly enhance these values. 
For example, within Minimal Supersymmetric extensions 
of the SM (MSSM), in the large tan/3 regime, &(BP S 
p + p~) receives contributions proportional to tan 6 /? ifTTl . 
where tan j3 is the ratio of the vacuum expectation values of 
the two neutral CP-even Higgs fields, and can differ signifi- 
cantly from the SM prediction. The LHCb analysis is done 
by clasifying B® — > p + p~ candidates in bins of a 2D pa- 
rameter space made by the invariant mass and a multivari- 
ate clasifier which condensates geometrical and kinemati- 
cal information of the event. The signal expectation in each 
bin is calculated using data from control channels such as 
59 v —> h + h'~ and B + —> J/i//K + . The background expec- 
tation is calculated by interpolating from mass sidebands. 
The 2?^ — > h + h'~ peaking background yield is calculated 
using 7i — > p and K —> p misidentification probabilities 
obtained from data using decays such as A — > pn~ and 
D° — > K + n~ . The signal and background expectations are 
compared with the distribution of observed events, and the 
limits are set using the CL S method H12I1311 . The !B(B° S -> 
p + p~) and £>(B° —> p + p~) upper limits obtained are: 

&(B° S -> p + p~) < 1.2 (1.4) x 10" 8 at 90 % (95 %) CL, 
£(B° -> /iV) < 2.6 (3.2) x 10" 9 at 90 % (95 %) CL. 

Fig. [2] shows the luminosity needed to impose stronger 
limits or to achieve a 3cr evidence of B° s — > p + p~. 



5 Conclusions 

As can be seen in Fig. [I] there is good agreement between 
recent SM predictions and LHCb's measurement of A FB , 
F L and dBF/dq 2 in the six q 2 bins. In a 1 < q 2 < 6 GeV 2 
bin, LHCb measures A FB = -0.10^+0.05, F L = 0.57^;}^ 
0.03 and dBF/dq 2 = 0.39 ± 0.06 ± 0.02, to be compared 
with theoretical predictions of A FB = -0.04^ ^, F L = 
0.74+°$ and dBF/dq 2 = (0.50+g;}J) x 10" 7 respectively. 
The experimental uncertainties are presently statistically 
dominated, and will improve with a larger data set. Such a 
data set would also enable LHCb to explore a wide range 
of new observables lfl4l . 

In 340 pb -1 of pp collisions at a centre of mass energy 
of - 7 TeV the most precise measurement of !B(B° — > 
cpy) has been performed, giving: 



&(B ^ K*°y) 
®(B° s ^<Py) 



1.52±0.14(stat)±0.10(syst)±0.12(/ s // rf ) 

(3) 

The S(B° S -> p + p~) and S(B° -> p + p~) upper limits 
obtained by LHCb are: 

&(B° S -> V) < 1.2 (1.4) x 10" 8 at 90 % (95 %) CL, 
£(B° -> p + p~) < 2.6 (3.2) x 10" 9 at 90 % (95 %) CL. 

In Fig. [2] the luminosity needed for a 3<x evidence as 
a function of tB(B° s — > p + p~) is shown. Approximately 
~ 2fb _1 are needed in the case that the value is equal to 
the SM prediction, but statistical fluctuations can make it 
possible with ~ 1 fb -1 . Fig.[2]also shows that exclusions of 
£>(B° S — > p + p~) down to the (2x) SM level would impose 
important constraints in region around the current NUHM1 
best fit point 031 . All the results presented here are current 
world best. 
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Fig. 1. A FB , F L and the differential branching fraction as a func- 
tion of q 2 in the six Belle q 2 bins. The theory predictions are 
described from Ref. Q. 
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Fig. 2. Luminosity needed in order to get afij ^ \i + \r 3cr 
evidence (top) or a 95 %CL exclusion in the presence of a SM 
signal (center). The bottom plot shows how upper limits in the 
10~ 9 level would constraint the region around the minimum of 
the NUHM1 fit from (B). 



